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It is found that in addition to the expected Kosterlitz-Thouless phase transition this 
model exhibits an infinite series of phase transitions at special values of the lattice spacing 
e pq = sm(np/2q). An unusual property of these transitions is that they are totally invisible 
in the double scaling limit. A method which allows us to explore the transition regions 
analytically and to determine certain critical exponents is developed. It is argued that 
phase transitions of this kind can be induced by the interaction of two-dimensional vortices 
with curvature defects of a fluctuating random lattice. 
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1. Introduction 

After a period of considerable progress in the study of two-dimensional quantum 
gravity and noncritical strings [|l| several fundamental problems in the field remained open. 

The most prominent among them is, perhaps, the problem of the c = 1 barrier @0. 
Indeed, string theories with the central charge of the matter c < 1 are by now rather well 
understood. On the contrary, very little is known about c > 1 strings, but their physical 
properties are probably quite different ||. 

The problematic nature of c > 1 theories can be partly attributed to the presence of a 
tachyon in their spectra. However, this is not the main obstacle in studying them. Indeed, 
in the random matrix model approach to string theory the cause of the difficulties appears 
to be more technical ||. In fact, it is perfectly possible to write down a matrix model 
which in the large iV limit reproduces the genus expansion of a .D-dimensional noncritical 
string with D > 1. What is lacking are the analytic tools necessary to explore such models. 

The reason is that, compared to c < 1 models, string theories with c larger than one 
have significantly more degrees of freedom. In the continuum formulation these additional 
degrees of freedom correspond to the transverse modes of the string, absent for c < 1. 
The appearance of new modes is especially clear in the matrix model language. The 
dynamical variables of matrix models are Hermitian N x N matrices, each characterized 
by A 2 independent matrix elements. If we diagonalize a Hermitian matrix, representing 
it as M = UXU*, these iV 2 degrees of freedom can be split into iV eigenvalues forming the 
diagonal matrix A and A 2 — A "angular variables" encoded in the unitary matrix U. In 
c < 1 matrix models the angular degrees of freedom decouple (or can be easily integrated 
out), leaving us with a theory of iV eigenvalues which contains all the information about 
the original string theory. For matrix models of D > 1 strings this no longer occurs. As a 
result, in more than one dimension the dynamics of the theory depends on all of A 2 ^> A 
variables. 

Mathematically, when the angular variables do not decouple the partition function of 
a matrix model involves a nontrivial integral over unitary matrices U. Explicit evaluation 
of such integrals would be useful for applications not only in matrix models but also in 
lattice gauge theories. In fact, the large A limit of QCD can be interpreted, through the 
Eguchi-Kawai reduction, as a certain integral over four unitary matrices . Thus even the 
large A QCD can be viewed as a very complicated multimatrix model, where the unitary 
degrees of freedom play a crucial role. 

In this paper we shall explore the effects of such degrees of freedom in the simplest 
model where these effects are nontrivial. The theory to be considered is a one-dimensional 



infinite chain of Hermitian matrices M n defined by the partition functional 

Z = I nrfM n exp{-iVTV g \ {Mn+ \~ Mn) \ eV(M n )\\. (1.1) 

n= — oo n= — oo ^ ^ ^ 



This model is, in fact, quite interesting by itself It represents a one-dimensional 

string theory whose target space, instead of being continuous, consists of an infinite number 
of equidistant discrete points. Furthermore, as we shall discuss below, the large iV limit of 
this model is related by a duality transformation to the two-dimensional 0(2) nonlinear 
sigma model coupled to quantum gravity or, equivalently, to one-dimensional bosonic string 
theory compactified on a circle of radius R = 1/e. 

The multimatrix chain exhibits a phase transition induced precisely by the unitary 
degrees of freedom discussed above Indeed, in the limit of infinitely small lattice 

spacing e —>■ our model describes the one-dimensional string theory with a continuous 
target space. Such a theory has c = 1. On the other hand, when the lattice spacing gets 
very large the interaction of neighboring matrices Tr(M n — M n+ i) 2 /2e becomes negligible 
compared to eTrV(M). Consequently, at e = oo the multimatrix chain decouples into a 
product of infinitely many identical one-matrix models. For cubic V(M) these one-matrix 
models represent pure two-dimensional quantum gravity, a system with c = 0. Therefore, 
the central charge of string theory corresponding to ( |1 . 1| ) changes from one to zero as e is 
increased. In fact, there is evidence that this change occurs sharply at a certain value of 
the lattice spacing e = e cr where the model undergoes a phase transition. 

Of course, being one-dimensional the matrix chain is substantially simpler than any 
of the prospective c > 1 models. However, we shall see that already in this simple theory 
the unitary variables produce new effects not encountered at c < 1. Surprisingly, we find 
that the phase transition at e = e cr is not the only critical point of ( |1 . 1| ) . It turns out 
that in addition this model has an infinity of critical points at e = e cr sm(Tcp/2q) labelled 
by positive integers p and q, where certain observables, such as the eigenvalue densities 
of matrices M n , develop singularities. These singularities are of universal nature (that is, 
they are the same for various potentials V(M) within a certain class) and evidence phase 
transitions. It is quite unusual, however, that none of these transitions affects the double 
scaling limit of the theory so that ( |1 . 1| ) has c = 1 for any e < e cr . 

It is rather amusing that the additional critical points fill densely the interval [0, e cr ]. 
This shows that once the unitary degrees of freedom are taken into account the system 



1 In this formula V(M) is a polynomial potential such as, for example, m 2 M 2 /2 + gM 3 /3 or 
m 2 M 2 /2 + c/M 4 /4. 



can become very complex. Perhaps, it is not entirely unlikely that a similar kind of phase 
structure might also arise in c > 1 matrix models. 

In addition to being on the border between c < 1 and c > 1 theories, there is another 
reason why the matrix chain is interesting. It models, in the large N limit, one-dimensional 
string theory compactified on a circle of radius R = 1/e. When defined on a discretized 
worldsheet this theory automatically contains vortices — topologically nontrivial configura- 
tions where the string winds around the target space circle as we follow the boundary of 
an elementary worldsheet plaquette. These vortices are suppressed thermodynamically for 
large R (or, equivalently, small e) but become favored as R is decreased. As a result, they 
induce a phase transition in very much the same way the vortices of the two-dimensional 
0(2) model drive the Berezinsky-Kosterlitz-Thouless phase transition [iQfl. In fact, this 
Kosterlitz-Thouless-type transition is precisely the transition at e = e cr which separates 
the c = 1 and c = phases of our theory. 

To explore the dynamics of the Kosterlitz-Thouless phase transition in string theory 



is a longstanding problem pjjP]||ll||. In this paper we shall present a computation which 



permits a systematic treatment of the transition region. In particular, we find that the 
eigenvalue density at the Kosterlitz-Thouless critical point has an unusual logarithmic 
singularity! p{x) oc \x\/ log(l/A|x|). However, we also find that the inclusion of lattice fluc- 
tuations gives rise to new Kosterlitz-Thouless-type transitions which would not occur on 
a regular flat lattice. These are the transitions at e = e cr sm.{iip/2q). They arise when the 
fluctuating random lattice develops a curvature defect somewhere around the vortex core. 
If such a defect introduces negative curvature, the energy of a vortex decreases. As a con- 
sequence, the effective temperature of the Kosterlitz-Thouless transition or, equivalently, 
the critical value of e for such vortices becomes lower than e cr . 

It is necessary to emphasize that these additional transitions do not modify the dou- 
ble scaling limit of the matrix chain. Therefore, they do not correspond to any new phase 
transitions in continuum one-dimensional string theory. However, they do exhibit univer- 
sality properties and this makes them worthy of consideration. Indeed, there are examples 
of phase transitions which appear as lattice artifacts in one model but play an important 
role in the continuum limit of another. For instance, the Gross- Witten phase transition 



in the one-plaquette model [12 1, which represents a lattice effect, is in the same universal- 



ity class with the Douglas-Kazakov phase transition [ 13 1 occuring in continuum large N 
Yang-Mills theory on a two-dimensional sphere. 



2 This should not be confused with the logarithmic singularity in the free energy of c = 1 
models where p(x) oc |x| without any logarithms. 



Since the analysis presented below is somewhat lengthy we shall first summarize the 
main steps of our computations and state the results. This shall be done in the next 
section. In sections 3 and 4 we shall present the details. In section 5 we check our results 
by studying the case of a special interaction potential V(M) for which the model can 
be solved exactly in terms of elementary functions. Finally, the vortex interpretation of 
defect- induced phase transitions shall be discussed in section 6. 



2. Large N Expansion of Infinite Matrix Chain 

The free energy of the infinite one-dimensional random matrix chain represents, in the 
large iV expansion, the string perturbation series for one-dimensional string theory with a 
discrete target space [U0. To see this consider a cubic V(M) = m 2 M 2 /2 + gM 3 /3 and 
expand ( |1 . 1| ) in powers of g. Such expansion is a sum of all Feynman diagrams with 
vertices and propagators which can be inferred directly from ( |1.1|) , 

(M n , af) M n ,„ s ) = N- l 5 Pl 5 aS D{n - ri), (2.1) 

where a,/?, 7 and 8 are matrix indices and the coordinate space propagator D(n — n') is 
given by 

7T 

D(n-n')= ( ^F. e i p( n - n ') - (2 2) 

V i 2tt e 2 m 2 + 4sin 2 (p/2)' ' 

— TT 

According to these rules, the contribution of any individual Feynman graph V with V 
vertices and G handles has the form 

+00 

n 1 ,...,n v = — oo (ij) 

As usual, the product of propagators goes over all links of the graph (ij) whereas the 
integers ni, . . . , ny refer to the vertex positions in the one-dimensional discrete coordinate 
space. In the string theory language these integers parametrize the embedding of the graph 
T into the discretised target space of the string. 

The weight assigned to a given graph in this theory, the product of propagators D(rii — 
rij), is not exactly equal to the discretized Polyakov string weight. Indeed, to obtain the 
Polyakov weight exp(— Sp[rii\) with 

SpN = ^e 2 K-%) 2 (2-4) 

(ij) 



one would have to choose the propagator Dp(n) = exp(— e 2 n 2 /2). Such a propagator could 
be generated only by a nonlocal, hard-to-deal-with matrix model. Fortunately, it is possible 
to argue @@ that the models with these two weights are in the same universality class. 
The reason is that both D and Dp have the same infrared behavior in momentum space, 
D(p) oc 1 — p 2 /p 2 . Since in one dimensional theory there are no ultraviolet divergences the 
replacement of Dp by D modifies only the short-distance, nonuniversal properties of the 
model. 

Alternatively, the large N limit of the infinite matrix chain ( p. . 1| ) provides a description 
of one-dimensional bosonic string theory compactified on a circle . However, in contrast 
with the discretized target space picture, this second interpretation is restricted to the 
leading order of large N expansion. That is to say, our matrix model reproduces correctly 
only the first term of the genus expansion for the compactified string. To explain why 
let us perform a duality transformation on the Feynman amplitude ( |2.3| ). As a first step, 
insert the momentum space representation ( |2.2| ) for each D in ( p.3| ) and do the sums over 
rrii using the identity 

+00 +00 



e ipn = 2n S(p + 2ttI). 



The result is 



z=— 00 



+00 w . , v 



g v N 2-2G £ Yl J DM J] [^(PL ~ 2nl k )] (2.5) 

l 1 ,...,l v = -oo (ij) _ 7r k=l 

where pij stands for the momentum flowing along the link (ij), D(p) is the momentum 
space propagator and p^ ot = J2j Pkj equals the sum of all momenta entering the vertex 
number k. Note that the vertices conserve momentum only modulo 2-7T. 

If r is topologically spherical it is easy to interpret ( |2.5|) in terms of compactified string 
theory. To this effect we replace the link momenta pij in ( |2.5|) by the string variables X a , 
to be defined shortly, which will be associated with the vertices of the dual graph T*. 
Since we consider a compactified string the variable X must live on a circle, so that X and 
X + 2nRn for any tiGZ are all identified. 

To construct X ai for each link (ij) G T we find a unique link (ab) of T* which crosses 
(ij) and impose the condition 

X a -X h = R PlJ . (2.6) 

If we fix the value of X at any one vertex (say, set X\ = 0) then equations (^]^) define 
X a as functions of p^ at all other vertices. In terms of these new variables ( |2.5| ) takes the 
form 



Fig. 1: A typical Feynman graph F generated by a matrix model (solid lines) and 
its dual r* (dashed lines). The discrete target space coordinate n is associated 
with the vertices of T while the compactified matter field X corresponds to the 
vertices of T* . 



which looks very similar to ( |2.3| ) except that the target space is now a continuous circle of 
radius R = 1/e. 

The equivalence we just described does not hold on higher genus graphs where in 
general ( [2.6| ) does not have a solution. However, for our purposes this does not matter. 
Indeed, below we shall restrict our attention to the large N limit of the chain model where 
only the spherical graphs survive and both interpretations apply. 

To explore the large N limit of the multimatrix chain we shall first integrate out 
the unitary degrees of freedom in ( |1 . 1|) and then use the saddle point method. It will be 
convenient to rescale all matrices M n = ^Ai n so that e disappears from the kinetic term 
in the partition function: 

/OO s OO N 

II dM n exp j ATtr £ [M n M n+1 - U(M n j\ \ (2.8) 

n= — oo n= — oo ^ 

where the potential U(M) is related to V(M), 

U(M) = M 2 + eV(y/eM). (2.9) 

As usual we shall diagonalize each matrix, Ai n = U n A n U^ with diagonal A n = 
diag(Ai jn . . . Xn,u) and express the Hermitian matrix measure dM. n in terms of U's and 
lambdas, 

dM n = A 2 (A n ) d\i y7l . . . d\ N , n dU n 



where A(A n ) is the Van der Monde determinant of eigenvalues A^ n and dU n refers to the 
Haar measure on the unitary group SU(N). Furthermore, it is useful to introduce the 
matrices V n = U^ +1 U n . 

Given this notation the partition function of the infinite matrix chain can be written 
in the form 



Z=j \[^ 2 {X n ) dX ltn . ..dX N , n dU n expjiVtr [v n A n V^A n+1 - *7(A n )] |. (2. 

n£2 n= — oo 



.10) 



The matrices V n represent the angular degrees of freedom in our model. Fortunately, since 
the one-dimensional lattice does not have closed loops, all V n are mutually independent 
and we can easily integrate them out. To do this one simply changes variables from {U n } 
to {V n } with the result 

Z= /n & 2 (K)dX 1 , n ...dX N>n dV n expjiVtr [v n A n V^A n+1 - U (A n )\\ . (2.11) 



n= — oc 



Note that no Jacobian arises when we pass from U to V. 

To integrate out V's one must first compute the following unitary group integral 

J dV exp [Ntv(VA n V^A n+1 )} = exp [iV 2 F(A n , A n+1 )] . (2.12) 



Then the resulting integral over the eigenvalues 
2 =J II A '( A -) dX ^ ■ ■ ■ dX N,ndV n exp jiV jr | iVF(A n , A n+1 ) - tr U(A r 

n£2 n= — cx 



(2.13) 



can be evaluated, at — >• oo, via the saddle point method. In fact, due to the translational 
invariance of the matrix chain the saddle point values of Aj, n are always independent of 
the site number n. The saddle point equations which determine these values are readily 
obtained by maximizing the integrand of (|2.13| ): 



2N 



<9F(A,A') 



dX % 



A'=A 



2 x ^ 



1 



N ^f. Xi - Xj 



0. 



(2.14) 



The unitary integral ( |2.12| ) has been computed in the classic paper of Itzykson and 
Zuber |T5|] . Unfortunately though, their result is not easy to use directly in the saddle point 
equation ( |2.14j ). Indeed, the Itzykson-Zuber formula expresses F through the determinant 
of a certain N x N matrix, a quantity which becomes quite complicated in the large N 
limit. 



As N goes to infinity it is convenient to characterize the eigenvalues composing the 
diagonal matrices A = diag(Ai, . . . , Ajy) in terms of the so-called eigenvalue densities p(A). 
By definition, Np(X) dX is the number of those eigenvalues among Ai, . . . , Xn which fall 
into an infinitesimal interval [A, A + dX] . To define the large N limit of ( |2.12| ) we enlarge 
diagonal matrices A n and A n _|_i in such a way that their respective densities of eigenvalues 
converge to the well defined smooth limits p n (X) and p n _|_i(A). Note that as a consequence 
of the above definition any density p(X) always obeys the constraint 

+00 

p(X)dX= 1. 



We shall now describe, without a derivation, how one computes F for infinitely large 
N. It turns out that the large N asymptotics of the Itzykson-Zuber integral ( |2.12|) is 
related to the classical mechanics of a certain one-dimensional integrable system |l6j . This 
system is a droplet of one-dimensional compressible fluid which has the following peculiar 
equation of state 



P 



-P 



(2.15) 



where both the local pressure P and the fluid density p depend, in general, on the one- 
dimensional coordinate x. The motion of such a droplet is described by Euler equations 



dp ^ d 

m + d-x [pv] 







dv dv 1 dP 
dt dx p dx 

= tt 2 — 
dx 



(2.16) 



where v(x) is the local velocity at point x. 

To compute F in equation ( |2.12| ) one seeks the solution of these equations which 
satisfies the boundary conditions 

p(x,t= 0) = p n (x) 
p(x,t= 1) = p n +l(x) 

where p n (x) and p n+ i(x) are the eigenvalue densities for the matrices A n and A n _|_i re- 
spectively. If the functions p(x,t) and v(x,t) yield such a solution then 

1 



(2.17) 



N 



+ 



lim F(A n ,A n+ i) 



r / <U 



dx p(x, t) 



v 2 (x,t) + ^-p 2 (x,t) 



^ J p n (x)x 2 dx- J J dx 1 dx2p n {x 1 )p n (x 2 ) \D.\Xl-X2\ 



(2.18) 



+ 



Pn+i(y) V 2 dy 



dy x dy 2 p n +i(yi)Pn+i(y2) In \yi - y 2 \ 



It is worth mentioning that this expression for F has a very special structure. Most 
importantly, the double integral 



S = — I <H I <lv /)(./:. t) 



v 2 (x,t) + ^-p 2 (x,t) 



(2.19) 



is actually the value of the classical action for the fluid droplet evolving according to 
equations ( |2TT6[ ) and ( |2T7| ). 

At this point we are ready to derive the large A" limit of the saddle point equation 
Q2.14]) . First of all, we need to know the derivative dF(A n , A n+ i)/d\i with respect to 
a given eigenvalue A^. At large A" such a derivative can be easily expressed (simply by 
manipulating the definitions) through the functional derivative of F[p n , p n +i] with respect 
to p n , 



N 



<9F(A n ,A n+1 ; 



d 5F[p n ,p n+1 ] 
dx Sp n (x) 



(2.20) 



x — Xi 



Furthermore, using ( |2.19|) and equations of motion ( |2.16| ) it is possible to check that 

d 5S[p n ,p n+1 ] 



and, therefore, 



dx Sp n (x 
d 5F[p n ,p n+1 ] 



-v(x,t = 0) 



2.21) 



v(x,t = 0) - x + 2 



p(y) dy 

x-y ' 



[2.22) 



dx 5p n (x) 

These relations are very natural. Indeed, the variation of the action S on any classical 
solution with respect to the initial coordinate of that solution (in this case p n (x)) always 
equals minus the corresponding canonical momentum. 

We can use ( |2.20|) and ( |2.22| ) to simplify the saddle point equation (|2.14| ) . Keeping in 
mind that for A^ — > oo 



f piv) dy 
J Aj 



V 



one easily reduces the saddle point equation to 

1 



v(x,t = 0) = —U'(x) — x. 
2 



(2.23) 



The mathematical setup of the problem is therefore as follows. Given the potential 
U(x) one is looking for the solution of Euler equations ( |2.16|) which satisfies the boundary 
conditions 

p(x, t = 0) = p(x, t = 1) 

. ( 2 - 24 ) 



v(x,t= 0) = -U'(x) 
2 



X. 



Once this solution is found the saddle point density p(x,t = 0) = p(x,t = 1) = p(x) can 
be used to compute Z or the free energy T = \ogZ. Indeed, for the quartic U(x) = 
m 2 x 2 /2 + gx 4 /4 it is easy to show that, as a consequence of ( |2.13| ) and ( |2.14| ) 

^ = -- p(x)x 4 dx (2.25) 



dg 4 

which immediately yields T(g) if p(x) is known. 

To say the same differently, imagine a fluid droplet of spatially dependent density 
p(x) (to be determined) being pushed with known spatially dependent initial velocity 
v(x) = U'(x)/2 — x. Demand that after one unit of time the density of that droplet evolves 
into the same p(x). This condition fixes p(x), hopefully uniquely, for any given v(x). 

Most remarkably, Euler equations (|2.16|) are explicitly integrable. This becomes clear 
if we combine v and p into one complex- valued function / = v + in p. One discovers that 



the two equations (|2.16|) are equivalent to the following complex equation on / 



commonly known as the Hopf equation. The Cauchy problem for the Hopf equation is 
easily solvable. If fo(x) — f(x,t = 0) is the initial value of the function / then the value 
of f(x, t) can be found from the implicit equation 

f(x,t) = f [x-tf(x,t)]. (2.27) 

Now we can impose the boundary conditions ( |2.24| ) in a more explicit form. By manip- 
ulating the implicit solution Q2.27 ) it is possible to show [16] that ( 2.24 ) gives rise to a 



functional equation on p(x). This functional equation is most easily formulated in terms 
of two auxiliary functions 

G±(x) = -U'(x) ± mp(x) (2.28) 
2 

and reads 

G+[G-{x)]=x. (2.29) 

It provides a direct relation between the interaction potential and the large N eigenvalue 
density of the infinite matrix chain. 



Unfortunately, little is known about solutions of such functional equations ||17|| . Cer- 
tainly, it seems difficult to solve ( |2.29| ) exactly for an arbitrary U(x). This, however, is not 
always necessary. Indeed, for most purposes one is interested only in the universal proper- 
ties of p -those independent of the particular U(x) chosen. For example, if a matrix model 



undergoes a phase transition the eigenvalue density p(x) usually develops a singularity at 
the transition point. Such singularities are typically of the form p(x) oc \x\ 5 (if the singular 
point is at x = 0) characterized by a universal exponent 5. In fact, 5 is closely related to 
the string susceptibility 7 str = 1 — 5, a quantity which plays an important role in string 
theory. 

Therefore, one might try to evaluate such exponents without finding the whole p. 
Luckily, equation ( |2.29| ) is well suited for just that. For concreteness and to make our 
problem well defined we shall concentrate on even potentials bounded from below, such as 

U(x) = -^f+ 9 -^. (2.30) 

Any other (say, cubic) potential can be treated the same way. To set the notation let us 
choose V(M) in the original definition (|1 . 1| ) asi 

V(M) = -2M 2 + gM 4 /4. (2.31) 

This, according to ( |2.9| ) means that 

2 = 2(2e 2 - 1) 

(2.32) 




For any given fixed e the multimatrix chain can be in one of the two phases, depending 
on the value of g. For large g (strong coupling) the height of the hump in the potential 
V(M) is relatively low so that the eigenvalue density p{x) has only a small dip around 
x = 0. For very small g, however, the hump becomes high and the "eigenvalue droplet" 
splits into two disconnected pieces. In between these two phases at a certain g = g cr (e) 
there is a phase transition where p(x) just vanishes at x — ► 0, 

p(x) oc |x| 5(e) . (2.33) 

The critical properties of the matrix chain encoded in 8(e) are not the same for various e. 
In fact, this model possesses a number of multicritical points — special values of epsilon— 
where 5 or other exponents can abruptly change values. One of the problems we shall 
address below is to find and investigate such multicritical points. 

Before we proceed further it would be convenient to summarize our findings: 



For a general quartic double well potential we can always rescale M and e in (1.1) to make 
the coefficient at M 2 equal —2. The physics of the model is not affected by such rescaling. 




X 



Fig. 2: A typical double-well potential V(x) and the corresponding eigenvalue 
density p(x) in the c = 1 matrix model at the critical coupling g = g cr (thick solid 
line), in the weak coupling phase g < g cr (dotted line) and for strong coupling 
g > g cr (dashed line). 

(a) If < e < 1 then 5(e) = 1, an old result established originally by Gross and 
Klebanov [[7]. This value of S is typical of c = 1 theories. That is to say, the discrete 
structure of the target space does not affect the string theory so long as the target space 
lattice spacing is less than one. 

(b) The Kosterlitz-Thouless point e = 1 is special. It separates two regimes where the 
infinite matrix chain has different values of 5 and describes different string theories. We 
find that for e = 1 the critical eigenvalue density (that is, the density at g = g CI ) develops 
a logarithmic singularity 



However, the universal properties of the matrix chain at the Kosterlitz-Thouless point 
are not exhausted by that. As we shall prove, an infinite number of corrections to this 
logarithmic singularity also happen to be universal. More precisely, the critical density 
p(x) can be expanded in the form 




(2.34) 




(2.35) 



q=0 p=l 



where y = log[l/(A|x|)] and I = logy. The numbers a pq which arise as coefficients of 
this expansion are entirely independent of the interaction potential V{M). The potential 
enters only through the nonuniversal scale A and through the terms of order 0{x 2 ). We 



shall also see that the series in (|2.35|) can be summed. The result is 

p(x) = \x\ a(x) (2.36) 
where cr(x) is determined by the transcendental equation 

- + loga = logf-^V (2.37) 



a \A\x , 

Quite remarkably, nothing in this equation except for A depends on the matrix model 
potential. Therefore, equation ( |2.37| ) ought to have a continuum, string-theoretic interpre- 
tation. In the continuum language the phase transition at e = 1 is induced by topologically 
nontrivial string configurations — vortices, which become strongly coupled for e > 1. From 
this point of view (|2.37 ) might contain information about vortex dynamics in the region 



where the "vortex gas" ceases to be dilute. 

(c) Surprisingly, the linear matrix chain has other special points. To observe them it 
is necessary to evaluate corrections to ( |2.33|) . For a generic e e]0, 1[ these corrections are 
of the form 

p(x) = |x|{ai(e) + a 2 (e)x 2 + a 3 (e)x 4 + . . . } (2.38) 

where the coefficients afc(e) are fixed uniquely by Euler equations. However, if e = 
sin(7rp/2g) = e pq we discover that the coefficient a q (e) goes to infinity. At these values 
of e the corrections to p(x) cannot be expanded in a power series. Indeed, we shall see that 
at e = e pq the correct expansion for p(x) contains a subleading logarithmic singularity 

p(x) = \x\ jai + a 2 x 2 + ... + aq-xx 2 ^-^ + a^x 2 ^" 1 ) logf^^J + . . .J. (2.39) 

These logarithmic singularities are universal. That is to say, they remain present if the 
quartic interaction is replaced by any other generic polynomial potential V(M). It is also 
possible to construct a number of special, "fine-tuned" potentials for which some of the 
singularities get eliminated. 

Mathematically, such logarithmic singularities are quite analogous to the Kosterlitz- 
Thouless singularity at e = 1. For example, if we shift g away from g CI (e) the shapes of 
p(x) at e = 1 and at e = sin(7rp/2g) are given by very similar expressions. The qualitative 
pictures of what happens at these points seem to be related as well. Indeed, there is 
evidence that the singularities at e = e pq , like the Kosterlitz-Thouless phase transition at 



e = 1, are due to effects of vortex proliferation. Such proliferation is, of course, impossible 
for any e < 1 in flat spacetime or on a regular lattice. On a random lattice the situation 
is different. There a pair consisting of a vortex and a negative curvature defect centered 
close to each other may become a favored configuration even when e < 1 . We shall see that 
although such pairs do not make a leading singular contribution to the eigenvalue density 
or the free energy of the model, they are still sufficient to cause a mild singularity in p(x). 

(d) If e > 1 and V(M) is quartic the Euler equations are consistent with 6(e) = 2. This 
is the same value of 5 one observes for e — > oo when the matrix chain becomes a collection 
of infinitely many decoupled one-matrix models. The transition between 6 = 1 at e < 1 to 
5 = 2 for e > 1 has an interesting interpretation in terms of the hydrodynamic picture. It 
turns out that for e > 1 the motion of the liquid droplet prescribed by ( |2.16|) and ( |2.24j ) 



results in the formation of a shock. That is, at a certain moment in time, t s h = l/2e 2 the 
spatial derivative of p(x, t) at x = becomes infinite. This phenomenon never happens for 
e < 1. Moreover, the droplet picture provides a natural order parameter for the Kosterlitz- 
Thouless transition. This order parameter is defined by £(e) = p(x = 0, t = 1/2) where 
p(x,t) is the solution corresponding to g = g CT (e). The quantity ((e) has the property 

C(e) = if e < 1 

(2.40) 

C(e) > if e > 1 

and will play an important role in our analysis of the e > 1 phase. 

(e) Finally, we shall study the character expansion of the infinite matrix chain. The 



idea is to expand the integrand of the partition function ( 2.11 ) in a certain sum over U(N) 



characters Then within each term of this sum the unitary degrees of freedom V n 

can be integrated out with ease. At large iV the resulting sum over representations of 
U(N) is dominated by a single "saddle point representation." The highest weights of this 
representation (m, . . . , njv) are all integers of order N. It is conventional to characterize 
such highest weights by the "Young tableau density" pi(h) defined as the density of points 
hi = (rii — i)/N + 1/2 in a small interval around h. Note that since the highest weights 
are ordered the Young tableau density can never be greater than one. We shall prove that 
the Young tableau density pi(h) is directly related to the density pi/2(x) = p(x,t = 1/2) 
which arises in the hydrodynamic picture: 



npi/2 



COS 



irpi(h) 



2vW'^Y (2.41) 



This relation shows that there is a lot in common between the Kosterlitz-Thouless 
phase transition and the Douglas-Kazakov-type transitions which are observed in two- 



dimensional QCD or in dually weighted graph models [13] [18]. In the Douglas-Kazakov 



transitions the Young tableau density is always less than one for weak coupling but devel- 
ops a plateau with pi(h) = 1 for any h G [—/&*, h*\ at strong coupling. Exactly the same 
happens in the matrix chain. Indeed, since for h = h* the density pi(h#) = 1, equation 
fl2.41|) predicts 

C = Pi/ 2 (0) = -v / ^- (2-42) 

Equation Q2.40Q now implies that the plateau in pi is present only for e > 1, in complete 
analogy with the Douglas-Kazakov transition. The physical pictures of the two transitions 
are also analogous. Indeed, both the Kosterlitz-Thouless and the Douglas-Kazakov tran- 
sitions are induced by the topologically nontrivial configurations (vortices and instantons, 
respectively) [|13|][[nj which dominate in the strong coupling phases of two-dimensional 
0(2) model and two-dimensional QCD. 

A more detailed analysis of the infinite matrix chain at e > 1 will be given separately 
in another report. Below we shall discuss the Kosterlitz-Thouless phase transition and the 
subleading logarithmic singularities which arise in the e < 1 phase of the model. 



3. Critical Behavior of Eigenvalue Densities for e < 1 

In this section we shall explore the critical properties of the eigenvalue density for 
lattice spacings smaller than one. 

Let us first concentrate on the case when the coupling constant g is exactly equal to 
the critical value <7 cr (e). In the language of continuum theory this would correspond to 
setting the renormalized cosmological constant to zero. Furthermore, since the singularity 
of p(x) occurs at x = we can restrict our attention to small values of x. 

For quartic U(x) given by (|2T30|) and ( ggg ) the boundary conditions (ggg) read 

v(x, t = 0) = vq(x) = -U'(x) — x 

Naively, when x is small we can neglect the cubic term in v and consider the simplified 
boundary problem where the initial velocity is given by 

v (x) = -2e 2 x. (3.2) 

The solution of Euler equations subject to this boundary condition is easy to construct. 
We simply observe that the ansatz 



irp(x, t) = a{t)\x\ 
v(x, t) = (3{t) x 



(3.3) 



is consistent with Euler equations (|2.16|) provided 



(3.4) 




where the dot denotes the time derivative. 

These equations can be easily solved. One introduces the complex valued function 
fi(t) = (3{t) + ia(t) which, as a consequence of ( |3.4f) obeys 

fi + fi = 0, (3.5) 

so that 

flit) = T^— (3-6) 

t — Co 

with a certain complex valued constant cq. This constant should be determined from the 
boundary conditions. The first condition in ( |2.24j ) amounts to lm/i(0) = Im/i(l) which 
entails 

1 

CO = - + IT 

with real tq. Then from the condition on the initial velocity, Re/i(0) = — 2e 2 one deter- 
mines the constant tq 

t = ±VT=? (3.7) 
and, finally, the eigenvalue density at small values of x 



np(x,t = 0) = a(0)|x| 



2eVl - e 2 Id. 



(3.f 



The critical exponent for such p is obviously 5(e) = 1. Furthermore, we see that the 
expression for p{x) given by (|3.8| ) degenerates when e = 1. This is the first sign of the 
Kosterlitz-Thouless phase transition which occurs at that point. 

The simple solution we just constructed exhibits several important general features. 
First, it is symmetric with respect to time reflection around t = 1/2: 

p(x, 1 — t) = p(x, t) 

(3.9) 

v(x, 1 — t) = —v(x, t). 



This property is in fact true for any solution of Euler equations which obeys boundary 
conditions ( 2.24 ). The reason is that both Euler equations and the boundary conditions 
are invariant with respect to the reversal of time. That is to say, if p(x, t) and v(x, t) yield 



a solution, so do p(x, t) = p(x, 1 — t) and v(x, t) = —v(x, 1 — t). Furthermore, due to the 
translational invariance of the matrix chain the final velocity at t = 1 is opposite to the 
initial velocity vo(x), 

v(x,t= 1) = -v(x,t = 0). (3.10) 

This, together with p(x,t = 0) = p(x, t = 1), implies that the solutions {p, v} and {p, v } 
coincide at t = 0. Consequently, these two solutions coincide also at all later times, and 
(p|) follows. 

The relation between the initial and final velocities, which we needed for the proof, is 
almost obvious. Indeed, the action S[p n , p n +i] is a symmetric functional of the densities 
p n and Pn+i- Therefore if p n = p n +i the functional variation of S with respect to p n equals 
the variation of S with respect to p n +i- On the other hand, the end velocity is related to 
the variation of the action by 

d ^[Pn,Pn + l]^ +v ^ t=1) {3U) 



dx 5p n+1 (x) 



with a plus sign. On comparison with equation ( |2.21|) we immediately deduce ( |3.10|) . 
An important consequence of the time reflection symmetry is 

v(x,t = l/2) = 0. (3.12) 

In other words, after one half unit of time the fluid completely stops moving. This leads 
to the following qualitative picture of the droplet evolution. At the initial moment of time 
t = the fluid in the vicinity of x = is moving inwards. The pressure P = — 7r 2 p 3 /3 acts 
to slow down and stop this motion. If the initial p(x) has been chosen properly then at 
t = 1/2 the velocity of the fluid vanishes simultaneously everywhere. At larger t the forces 
of pressure make the droplet move again, repeating backwards in time the evolution from 
t = to t = 1/2. 

Let us now see how the functional equation ( |2.29|) reproduces the result of (|3.8| ) . This 
equation involves the density of eigenvalues taken at a complex point, p[G-(x)}. Such an 
object is certainly well defined if p(x) is analytic. However, the density profile of interest 
to us, p(x) oc \x\ cannot be continued analytically into the complex plane. 

The resolution of the arising difficulty is quite simple. One should throw away the 
absolute value sign and consider instead p(x) oc x. Such prescription does yield the correct 
answer. Indeed, if ixp(x) = ax the functions G + and G_ are given by 



G±(x) = (l-2e 2 ±ia)x. 



(3.13) 



Then equation ( |2.29|) reduces to 



(1 - 2e 2 ) 2 + a 2 = 1 
which predicts the correct value of a 

a = 2e\/l-e 2 . 

This prescription works because no fluid ever flows through x = 0. In fact, the regions 
x > and x < do not really interact with each other. Furthermore, nothing in Euler 
equations requires the density to be manifestly positive. For instance, we can flip the sign 
of p(x,t) for x < 0. The density so obtained, p(x,t) oc x will be analytic in x and will 
still satisfy Euler equations. On the other hand, the x > part of the eigenvalue density 
remains unchanged allowing us to read off the correct answer. 

So far the quartic interaction gtrM 4 /4 has not played any role in our analysis. An 
identical effect occurs in the standard formalism of c = 1 string theory based on matrix 
quantum mechanics @[]7J]. There, too, one expands the potential around the local maxi- 
mum and only the quadratic terms matter. In the hydrodynamic picture the interaction 
terms modify the initial velocity of the droplet by an amount negligible for small x. There- 
fore, for x — > the effects of interaction can be taken into account via perturbation theory. 
The validity of such perturbative expansion is ensured by the smallness of x even though 
the coupling g is of order one. In a sense, we are lucky that this perturbation theory is 
valid precisely in the critical domain. 

When the lattice spacing gets bigger than one such a simple picture ceases to be valid. 
For e > 1 imparting the initial velocity v = — 2e 2 x to a droplet configuration p{x) oc \x\ 
will always make the droplet collapse. In this case the cubic correction to the initial v 
given by gx 3 /2 provides a small, but essential, amount of outward directed velocity that 
is needed to prevent such collapse. 

Quite surprisingly, the interaction effects are not completely trivial even if e is less 
than one. This can be seen already in the first order of perturbation theory. Indeed, let 
us evaluate the first correction to the linear ansatz ( |3.3|) . To this end we expand p and v 
in powers of x 

7rp(x, t) = ai(t)\x\ + a2(t)\x\ 3 + . . . 

(3.14) 

v(x,t) = fait) x + /3 2 (t) x 3 + ... 

Technically, it is more convenient to work directly with the Hopf function / = v + iirp. For 
x > this function has a power series expansion in x, 



f(x,t) = f 1 (t)x + h(t) x 3 + 



(3.15) 



Substituting this into the Hopf equation (|2.26| ) and collecting the terms of order x 3 we 
obtain an ordinary differential equation on f2(t) 



/ 2 + 4M = 0. (3.16) 



Now we can use the explicit solution for fi given by (|3.6|) to find 



f2(t) = 102 (3.17) 
(r - zr ) 4 

where r = t — 1/2 and C2 is a constant of integration. At r = we must have v = which 
means that C2 is real. Further, the boundary condition ( |3.1|) yields Re/2^ = 0) = g/2 
fixing C2 to be 

9 9 / Q -| Q\ 

° 2 ~ ~32e 4 sin(4arcsine) ~ 128e 5 (2e 2 - 1)VT^' ' 
As a result, we find the first correction to the eigenvalue density 



g\x 



3 



7T P (x) = 26 y/l - e 2 |x| - _ + g ' . + • . . . (3.19) 

2 tan(4 arcsme) 

If e = l/v2 this correction becomes infinite. The technical reason for that is very clear. 
When e = l/y/2 the value of /2(t) at t = is purely imaginary, /2(C)) = — 4ic2- Therefore, 
even though we do perturb the eigenvalue density by x 3 terms, no change in the initial 
velocity arises. In this way, the x 3 mode at e = l/y/2 could be viewed as a "resonance" of 
the droplet. Certainly, a similar picture of resonances should also appear in the equivalent 
approach to c = 1 strings based on matrix quantum mechanics. 

Whenever such resonances occur the perturbative expansion needs to be modified. It 
turns out that an adequate modification is given by the following formula 

f(x, t) = h(t) x + f 2 (t) x 3 logx + j 2 {t) x 3 + . . . . (3.20) 

For this to be consistent with the Hopf equation f%{t) must satisfy the old differential 
equation ( ^.16| ) while f2(t) should obey 



/ 2 +4/ 1 / 2 + /i/ 2 = 0. (3.21) 

Consequently, the function f2(t) is still given by formula ( |3.17| ) with a real C2 and To = 
y/l — e 2 /2e =1/2. However, the constant C2 is not fixed by the boundary conditions on 
f2(t) anymore. Instead, it will be determined from the boundary conditions on f 2 . 



Since the initial velocity Vq(x) contains no logarithms we must demand that Re/2^ = 
0) = 0. According to ( |3.17| ) this is possible only if e = l/y/2. In other words, it is not 
always possible to add x s log x to the eigenvalue density — for that the lattice spacing must 
take a special value. And, in agreement with the picture of resonances, it is precisely at 
this special value when the original power series ( |3.19|) fails. 

Given an explicit solution for f2 (with C2 undetermined) we can find f2 from ( |3.21| ). 
Imposing the usual conditions Rej 1 ^ = 1/2) = and Rej 1 ^ = 0) = g/2 fixes C2 in ( |3.17| ) 
to be C2 = g/2iT and leads to the following expression for /2(£) 



hit) 



(3.22) 



with a yet undetermined constant c<i- The eigenvalue density which corresponds to such / 
has a logarithmic singularity at x = 



Tip(x) = \x\ + — |x| J log — — + C(x 4 ). (3.23) 

TV \X\x\J 

The scale A is related to C2 by C2 = (g/4n) log(A 2 /2) and cannot be determined from 
perturbation theory. To find such a scale, as well as the specific numerical value of the 
critical coupling g = <7 C r(e) one would have to construct the eigenvalue density for all, 
not only small, x and impose certain conditions on the global analytic structure and 
normalization of p(x). These conditions shall be discussed in more detail at the end of the 
next section. However, notwithstanding any ambiguity in the value of ci the structure of 
the singularity in p(x) is still fixed uniquely by the Hopf equation. 

Similar logarithmic singularities arise also in higher orders. In general, for e = e pq = 
sm(%p/2q) the eigenvalue density has a logarithmic correction of order |x| 2<3_1 log[l/(A|x|)]. 
This can be proved in a rather elementary fashion with the aid of the functional equation 
( |2.29|) . Indeed, if the eigenvalue density does not have logarithms the functions G±(x) can 
be expanded in a power series 

oo 

G±{x)=Y,bix 2k - 1 . (3.24) 

fc=i 

By construction, G+ and G- are complex conjugate, so that = (b^)*. Furthermore, 
since KeG±(x) = U'(x)/2 is a cubic polynomial all b^ with k > 3 are purely imaginary 
and thus satisfy b^ = —b^. 

Let us try to solve the functional equation ( |2.29| ) order by order in x. In the first 
order, bfb^ = 1 which means bf = exp(±i(p) with a certain real (p. Actually, comparing 



( |3. 24 ) to ( j3.13|) we see that the angle ip is related to the lattice spacing e by cos cp = 1 — 2e z 



or, equivalently, e = sin(</;/2). 

If we know all with k < q — 1 and wish to find b^ we expand (a;)] in powers 

of x and collect the terms of order x 2q ~ l . As a consequence of ( |2.29|) the sum of such terms 
must vanish. This yields 

W"' + bib- + ( deP ; ± nd } = 0. (3.25) 

9 9 [ only on b x , . . . , 6^ J 



For certain values of e equation (|3.25|) degenerates. Indeed, since b q = —6+ and b x 
exp(±ip) we find that 

6+(6r) 29_1 + b+b~ = -2ib+ sin(qp) e" 1 ^" 1 ^. 



If </? = 7175/(7 the variable 6+ disappears from the right hand side of ( 3.25 ) . Then, unless the 



rest of ( j3.25|) vanishes at the same time, we encounter a contradiction. As a consequence, 
for these values of (p or, equivalently, for e = sm(irp/2q) the eigenvalue density develops 
an additional singularity. We have checked explicitly that this is indeed true in the first 
seven orders of perturbation theory. 

The same type of reasoning demonstrates that the position, the order and the log- 
arithmic nature of arising singularities are preserved if the quartic interaction V(M) is 
replaced by a more complicated potential. There may be, however, exceptional situations. 
They occur when the terms enclosed by the curly brackets in equation ( |3.25 ) vanish si- 



multaneously with sm(qtf>). For that to happen the coefficients of the matrix potential 
V(M) must be adjusted in a special way. The eigenvalue density which corresponds to 
such fine-tuned, "multicritical" potentials will exhibit only some, or none, of the above 
singularities. For instance, in section 5 we shall construct an example of a nonpolynomial 
potential which gives rise to only one singularity at e = 1/V%- The model with such a 
potential can be solved exactly. In fact, it is a generalization of the Penner-type model 
where an exact solution can be gotten independently using the method of loop equations 

The above considerations apply when the coupling constant g strictly equals g C v(^)- 
Only then p(x) oc \x\ around x = 0. Whenever g is shifted away from <? C r(e) (or, in the 
language of continuum theory, when the cosmological constant is given a nonzero value) 
the eigenvalue density gets perturbed. Such perturbed shape of p(x) would be useful to 
know for several reasons. First, if one wishes to explore the critical behavior of the free 
energy, say using equation (|2.25| ), one needs to know p{x) in a certain interval of couplings 



around <? C r( e )- Further, we shall see that such analysis provides an independent derivation 
of the logarithmic singularities which were simply guessed in equation ( |3.20|) . 

For a small nonzero 5g = g~g C r{^) the |x|-type singularity in p(x) disappears. Instead 
the eigenvalue density develops a pair of close square root branch points. These branch 
points merge when 5g vanishes producing p(x) oc \x\. In fact, the analytic structure of 
exactly the same kind occurs in the c = 1 matrix model with a continuous target space 



where the eigenvalue density is given by the WKB formula H[21 



7ip(x) = \J2\_E — U(x)] (3.26) 

with U(x) = -k 2 x 2 /2 + gx 4 /4. 

Remarkably, the existence of two square root branch points is fully consistent with 
Euler equations. To see this let us consider the following "hyperbolic" ansatz for p and v 



np(x,t) = Jb{t) + a 2 it) x 2 

(3.27) 

v(x, t) = (3{t) x 

This ansatz is consistent with Euler equations provided that a(t) and (3(t) obey the old 
equations (|3.4|) while b(t) satisfies 

6 + 2/36 = 0. (3.28) 



Actually, the linear ansatz ( |3.3| ) that we considered previously is simply a 6 = case of 
the hyperbolic solution. 

The eigenvalue density which correspond to (|3.27| ) 



irp(x) = ^b(0)+4e 2 (l-e 2 )x 2 (3.29) 

has two square root branch points. Moreover, such p(x) can be viewed as a generalization 
of the WKB formula ( |3.26| ) to finite nonzero lattice spacings. Indeed, for small x the WKB 
solution has the form 



7rp(x) = v 2E + k 2 x 2 

identical to ( |3.29| ). The parameter b = 6(0) in the hyperbolic solution and the Fermi energy 
E in ( |3.26|) play extremely similar roles. Both 6 and E are in fact functions of 8g which are 
determined by the normalization condition that the total integral of the full, exact p(x) 
equals one. Obviously, for small coupling deviations Sg the parameters 6 and E are also 
small. 

In the critical region of interest to us both 6 and the eigenvalue magnitude x are small 
but can be of the same order with respect to each other. That is, if the characteristic 



width of the double well potential V(M) equals a then x ~ Vb <C a. It is exactly this 
region of eigenvalues — a small interval of size ~ vb around the top of the potential — 
that is important, for example, in the matrix quantum mechanics representation of c = 1 
strings. For such x and b equation ( |3.29|) yields a good approximation to p(x) which can 
be systematically improved upon. 

The corrections to ( |3.29|) can be easily computed using the functional equation 
G+[G-(x)] = x. We shall look for the corrected eigenvalue density in the form 



7rp(x) = y b + x 2 sin ip + gr(x) (3.30) 

where (p is defined by e = sin(<£>/2) and r(x) represents the first correction. For x 2 ~ b 
the quantity r(x) should be small of order b 2 or, equivalently, x 4 . Actually, for the special 
case of b = the estimate r(x) ~ x 4 follows already from (|3.19| ). The functions G + and 
G- are then 

gx 3 I 

G±(x) = x cosip H — - — \-iJb + x 2 sin 2 (p + gr(x). (3.31) 

The approximation we are constructing is formally the same as the power series expansion 
in g. Indeed, in such an expansion each extra power of g comes along with either b or x 2 . 
For example, gr(x) ~ gb 2 is much smaller than b + x 2 sin 2 (p which is of order b. The same 
applies to gx 3 /2 ~ gbx which is much smaller than x cos <p. In short, the true parameter 
of our expansion is gb ~ gx 2 C 1 so that we can formally expand in g even though g by 
itself is not small. 

Doing this in ( |3.31| ) to first order in g yields 

G ± {x)=G { l\x)+gH ± {x) 

with 

,(o) 



G± (x) = x cos ip ± i y b + x 2 sin ip (3.32) 
and 

H ± {x) = X -±i m . (3.33) 

z 2\/b + x 2 siii cp 



Then the functional equation (|2.29| ) expanded to 0(g 2 ) reduces to 



G'l 0) [GL 0) (x)] H- (x) + H+ [GL 0) (*)] = 0. (3.34) 

We now substitute for G±^ and H± the explicit expressions of ( |3.32| ) and ( |3.33| ). This 
produces the following linear functional equation on r(x) 



r(y) — r(x) = xy 3 — x s y + (x 4 — y 4 ) cos p> 



(3.35) 



where the letter y stands for the combination 



y = x cos f — iy b + x 2 sin ip. 

When x 2 ~ y 2 ~ b both sides of this equation are of order b 2 . It is therefore natural to 
look for a solution of ( |3.35[) in terms of a polynomial "homogeneous" in x 2 and b 

r(x) = Rtx 2 b + R 2 x 4 . (3.36) 

Upon a certain amount of computation one finds that this form of r is indeed consistent 
with the functional equation (|3.35| ) provided R\ and R 2 equal 

Rl 1 



2cos<^ 2(1 -2e 2 )' 

1 i ( 3 - 37 ) 

#2 = cos = — - 1 + 2e 2 . 

2cos<^ ^ 2(1 -2e 2 ) 

Expressions ( |3.36| ) and ( |3.37| ) generalize the perturbative expansion ( |3.19j ) to the case of 

a nonzero cosmo logical constant. Not surprisingly, these expressions also develop poles at 

e = l/y/2, where the polynomial form for r(x) ceases to be a solution of ( |3.35|) . However, 

even for e = l/y/2 the equation on r(x) can still be solved. Remarkably, the solution will 

produce a logarithmic contribution in r automatically without adding it by hand as we 

did before in (|3~20| ). 

When e = l/y/2 our linear functional equation becomes 

r(y) — r(x) = y 3 x — x 3 y 

(3.38) 

y = — % v b + x 2 

We are interested in the solution of this equation which has two square root branch points 
at x = ±iVb and, on the real axis, is an even function of x regular at x = 0. 

Such a solution is defined on a complex plane with two cuts as shown in the figure 
and can be found using the method of dispersion relations!!. Indeed, the imaginary part 
of r(x) on the edges of the cuts follows directly from ( |3.38|) and equals 

Im r(-iu + 0) = u(2u 2 - b) V V - b. (3.39) 



To see this, consider equation ( |3.38|) with a positive real x. Then y = —iu = —iyb + x 2 is 
positioned on the right edge of the lower cut. Furthermore, Imr(i) =0. As a consequence, 
Imr(y) = Im(y 3 x — x 3 y) leading immediately to (|3.39|) . 

Once the imaginary part of r on any of the four cut edges has been found the imaginary 
part on the other three edges can be determined by analytic continuation. Then we can 
write down a dispersion relation for the function r(x) to fix its values on the real axis. 



4 A shorter way to solve this equation, suggested by J. Goldstone, is to use the substitution x = 
\/b svahij). Nonetheless, the solution based on dispersion relations will be necessary to illustrate 
several important points later, in section 4. 
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Fig. 3: The cut structure of the eigenvalue density in the critical region and the 
contour C used in the dispersion relation. The same cuts and integration contour 
will arise in our analysis of the Kosterlitz-Thouless phase transition in section 4. 



To this effect consider the integral 

taken along the contour C shown in the figure. The large radius loops closing C at 
infinity do not make any contribution to this integral. This follows from the nature of our 
approximation. Indeed, let us choose the loop radius R to be much larger than b but still 
much smaller than the characteristic size of the potential V(M). For \z\ < R the function 
r(z) depends only on z and b and has the general scaling form r(z) = b 2 s(z 2 /b). However, 
for z 2 ^> b the leading asymptotic behavior of r(z) should be independent of 6, simply 
because r(z) has a finite limit when 6^0, given by (|3^3l) . Thus for large £ = z 2 jb the 
function s behaves as s(£) ~ £ 2 so that r(z) ~ z 4 , perhaps with logarithmic corrections. 
As a consequence for large z the integrand in (|3.40|) decays as 1/z 2 and the contribution 
from infinity equals zero. 

Therefore, the only parts of the contour that contribute to I(x) are the four cut edges. 
In fact, due to the square root nature of the branch points the real part of r on the cuts 



cancels out of the integral so that 
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du 
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(3.41) 



where 



8r(iu) = —5r(—iu) = —2iu(2u 2 — 6)V" 



is the jump of the function r(z) across the cut. 

Alternatively, I(x) can be evaluated by taking residues at z = and z = x. This 
yields 

I( X ) = 2ni r M- 2 J% + % + (3.42) 

X° \X° X^ X J 

the numbers fo,fi and r2 being defined by r(x) = fg + f\x 2 + r2X 4 + 0(x 6 ). This allows 
us to compute r(x). Evaluating the integral in (|3.41| ) and comparing the result to ( |3.42j ) 
we easily deduce 



r(x) = r + r\x 2 + r 2 x 4 (2x 2 + b) a/x 2 + b log 



y/x 2 + b + X 



X 2 + b — X 



(3.43) 



The new constants ro,n and r 2 which appear in ( |3.43| ) are related to the residues fo,fi 
and r 2 according to the formulas 



^0 =^0 

n =fi + 2b/ n 

r2 =f2 + 14/37T. 

At this point it is necessary to check that the expression we obtained satisfies the full 



functional equation (|3.38|) . The reason is that in the above analysis we threw out the real 
part of the functional equation taking only its imaginary part. It turns out that the full 
functional equation is obeyed if we impose an additional constraint r\ = r 2 b. Thus we 
finally get 



•(x) = tq + r 2 x 2 (x 2 + b) (2x 2 + b) a/x 2 + b log 



Vx 2 + b + 



X 



X 



(3.44) 



The parameter ro = r(0) can be set to zero by adequately redefining b. The parameter r 2 
plays the role analogous to C2 in fl3.23| ) and for the same reason cannot be determined by 
expansion methods. Note though that this parameter controls a totally regular, analytic 
contribution. The nonanalytic part of r(x) is, on the contrary, fully unambiguous. 

If we take b to zero we recover the logarithmic singularity in p(x) already found in 
( p.23|) . Indeed, given tq = the small b limit of (|3.44j ) becomes 

4 2x\ {Ax 2 \ 

r(x) = r Q + r 2 x log — — 

n \ b J 

which is fully consistent with ( |3.23| ) provided r 2 = — (2/ir) log(A 2 6/4). 

The expressions for the eigenvalue density we found in this section can be used, in 
principle, to determine the critical behavior of the free energy. We shall pursue this 
problem elsewhere. Furthermore, our expressions indicate that the analytic structure of 
the full exact solution for the matrix chain is rather complex. When viewed as an analytic 
function of the complex variable e this solution would have singularities at all points of the 
form e pq = sin(7rp/2g) (and, as a consequence, at all other e G [0, 1].) Remarkably, certain 
elliptic functions, like the Dedekind function 

oo 

V-\q)=q- 1 ^\{{l- q -)- 1 

n=l 

with q = exp(4zarcsine) have an analogous singularity structure. It would be quite inter- 
esting indeed if the exact p(x) for the matrix chain was actually connected to such elliptic 
functions. 



4. Kosterlitz— Thouless Phase Transition in Infinite Matrix Chain 

Mathematically, the Kosterlitz-Thouless phase transition in the matrix chain is quite 
similar to the logarithmic singularities we have discussed in the previous section. This 
transition occurs at the point e = 1 where the small x expansion for the eigenvalue density 
given by ( |3.8|) or ( |3.19|) ceases to make sense. 

In this section we shall demonstrate that the critical behavior of the eigenvalue density 
at e = 1 is given by a logarithmic law 

nix) oc t- (4.1) 

PK 1 log[l/(AM)]. 1 1 

Such behavior is somewhat atypical. Indeed, in matrix models the leading singularity in 
p(x) is usually powerlike, p(x) oc \x\ s . However, it is possible to prove that for e = 1 



a power law would be inconsistent with the saddle point equation. Although the explicit 
proof shall not be given here one can easily reconstruct it using the methods of this section. 

We can acquire some idea of what p(x) should be by inspecting the e < 1 solution 
( P.8|) . As e — > 1 the slope a in p(x) oc a\x\ tends to zero. This suggests that at e = 1 and 
x — > the density p(x) should vanish faster than a\x\ with any a. On the other hand, 
when e = 1 the asymptotics p(x) oc \x\ s with 5 > 1 contradicts Euler equations. Therefore, 
it is natural to look for p(x) which vanishes faster than \x\ but slower than \x\ 5 for any 
6 > 1. 

Expression has exactly this property and, most importantly, is consistent with 
Euler equations. The most convenient way to see this is to use the functional equation 
( |2.29| ). At e = 1 the functions G + (x) and G-(x) are 

3 

G±(x) = -x + — ±inp(x). (4.2) 

Such form for G± is very special. Indeed, if we neglect both the gx 3 /2 and mp{x) terms 
altogether the resulting functions G±(x) = —x would nonetheless satisfy G + [G-(x)} = x. 
From this point of view, the Kosterlitz-Thouless phase transition is just the most senior 
of the logarithmic singularities found in section 3. At these singularities, which occur for 
e = sin(7rp/2g), the small x behavior of G+ and G- is given by ( |3.13|) , or 



G±(x) = e ±l7Tp/q x + 



x3l ° g, i 



(4.3) 



If we raise G+(x) or G-(x) to the functional power 2q — that is, compose G+(x) with itself 
2q times, we shall get 



G+ G+[...G+(x)] ... 



x + . . . . 



2q times 

One could say that, up to the higher order terms, the functions G± are the fixed points 
of order 2q with respect to the operation of functional composition. All the logarithmic 
singularities of the matrix chain are in one-to-one correspondence with such fixed points. 
In this language, the Kosterlitz-Thouless phase transition is represented by the fixed point 
of the lowest — second — order. Consequently, there should be no surprise that at e = 1 we 
also encounter the logarithmic type of scaling behavior. 

Let us now add Tvp(x) = A\x\/ log[l/(A|x|)] with an unknown coefficient A and expand 
G+[G- (x)] around x = 0. In doing this expansion we shall use again the prescription which 
directs us to throw away all absolute value signs. The justification for that will be given 
later. In short, we write 

„ / x iAx 

G±{x) = -XT ; — tt-^ + • . . . 

log (Ax) 



The terms denoted by dots include the gx 3 /2 term coming from the interaction potential 
and the higher order corrections to the eigenvalue density. Generally, we expect these 
corrections to be down by powers of 1/| log(Ax)| or log | log(Ax)|/| log(Ax)|, like the higher 
order terms in the series (|2.35|) . Expanding the right hand side of the functional equation 
around x = yields 

iAG-{x) 



G+[G-(x)] = -G-{x) 



[x 



+ ... 



We can now replace log 
to get 

G+[G 



log [AG 

x — iAx log (Ax) + iAx log(— Ax) 
-Ax) = log(Ax) — Z7r and expand ( 
A(A 



iA 



x 
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log Ax 



1 (4-4) 

log(Ax) ~'~ 
in powers of 1/| log(Ax)| 



(4.5) 



log (Ax) 

Therefore, the functional equation G+[G-(x)] = x is indeed satisfied to the leading order 
provided A = tt, in complete agreement with ( |4.1| ). 

Note that in the last line of (|4.4|) we kept the terms of order x/ log 2 (Ax). At first 
sight this may seem redundant and even incorrect. After all, exactly such terms have 
been completely neglected in the original expression for p(x) given by (5T). Nonetheless, 
these additional terms change only the next, 0[x/\ log 3 (Ax) |] order of expansion. They 
cancel out from the xj log 2 (Ax) part of the series very much like the leading term itself, 
p(x) oc x/| log (Ax) | contributes only xj log 2 (Ax) and not x/| log(Ax)| to the right hand side 
of®. 

At this point we would like to show how one constructs the solution of Euler equations, 
in terms of v and p or the Hopf function /, which corresponds to ( f4.1| ). There are three 
reasons why this is required. First, it will provide us with a method to compute all 
logarithmic corrections to (|4.1| ), prove their universality and sum the correction series in 
( |2.35|) all at once. Second, we shall be able to justify our way of dealing with absolute 
value signs in (fj~4l). Third, functional equations like (|2.29|) can have spurious solutions 
which must be detected and eliminated. An example of such spurious solution is given by 
a power series 
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(4.6) 



with real coefficients a^. By expanding G+[G-(x)} in powers of x and equating the result 
to x one can find all in a recursive fashion without encountering apparent contradictions. 
Nonetheless, there is no solution of Euler equations which would give 



7rp(x, t = 0) = 7rp(x, t = 1) = 



a k x 



2k 



k=i 



while obeying the boundary condition v (x, t = 0) = — 2x+gx 3 /2. In other words, a solution 
of the functional equation (|2.29|) does not always yield a solution of the full hydrodynamic 
boundary problem ( |2.16| ),( |2.24j ). The reason is, when viewed as an analytic function of the 
complex variable x the eigenvalue density p(x) may have several different regular branches. 
One of them takes real values for real x and represents the actual eigenvalue density of 
the infinite matrix chain. The power series expansions we construct refer, by definition, 



only to this branch. However, the functional equation (|2.29|) involves the eigenvalue density 
computed at a complex point, p[G-{x)\ which might take us to an entirely different branch 
of p. The relationship between this other branch and the "proper" branch expandable in 
a series can in reality be very complicated. 

We can check whether a given solution of ( |2.29| ) suffers from branch choice problems 
by computing v and p at t = 1/2. A solution which is consistent with the complete set of 
hydrodynamic equations (|2.16|) , (|2.24|) shall have v(x, t = 1/2) = as required by the time 
reflection symmetry ( |3.12| ). Let us prove that this is indeed the case for the logarithmic 
eigenvalue density Q4.1|) . To this end we shall use the implicit solution of the Hopf equation 
given by formula ( |2.27| ). For e = 1 the initial velocity equals vq(x) = —2x + C(x 3 ) so that 
the initial value of the Hopf function is 



f (x) = -2x + 



17VX 



+ • 



log [l/(Ax)j 

Then the value of / at t = 1/2, denoted /i/2(x) can be determined from the equation 



(4.7) 
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x- 2/1/2(3) 



We are interested in solving this transcendental equation only in the limit of small x. It is 
easy to check that the relevant solution is given by 



4i ( 1 
f 1/2(2) = — xlogf j- ) + 



(4.9) 



where the dots stand for the subleading corrections of order 0[aTog | log(Ax)|]. Indeed, for 
such /1/2 
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where we neglected the real valued terms which involve double logarithms as well as real 
constant terms. The imaginary constants, on the other hand, can be kept for the reason 
explained after equation (f4.5[ ). In this approximation the left hand side of (|4.8[ ) equals 



2ix ( iri\ 2ix ( 1 

log x - — = x log - 

7T \ 2 J n \x 

which up to the terms of order C[xlog | log(Ax)|] coincides with x — /i/ 2 (x)/2. 

The expression for fi/ 2 in (|4.9|) is purely imaginary. Therefore, the fluid velocity at 
t = 1/2, given by the real part of /1/2, does vanish. Let us now examine the corresponding 
fluid density, 

4x / 1 \ 

7rpi/ 2 (x) = — logf j-\ + .... (4.10) 

Naively, we could expect that this expression receives logarithmic corrections, generally of 
order (9[a;log| log(Ax)|]. However, in a most remarkable way all such corrections happen 
to cancel. To avoid any misunderstanding, it would be false to say that the matrix chain 
eigenvalue density itself does not get logarithmic corrections. It does. But the quantity 
Pi/2{ x ) is still very simple, much simpler than p{x). In particular, the infinity of logarithmic 
terms in the double series (|2.35|) are fully summarized by equation ( |4.10|) . That is to say, 
the only corrections to Q4.10| ) are the 0(x 3 ) terms due to the quartic interaction in the 
matrix potential. 

To prove this we shall turn our arguments around. Imagine starting at t = 1/2 with 
a droplet of the density given by ( |4.10|) and zero velocity. Let the droplet move according 
to Euler equations. Then compute the density and velocity of the droplet at t = 1. 

Certainly, by time reflection we can always reverse the droplet motion to reconstruct 
the < t < 1/2 part of the trajectory. As a result the boundary condition p(x,t = 0) = 
p(x,t = 1) will be satisfied automatically. The only remaining boundary condition one 
will have to check is the equation on velocity 

3 

v (x,t = 0) = -v(x,t = 1) = -2x+ (4.11) 

2 



We shall now prove that, neglecting 0(x 3 ) terms, the logarithmic density of ( |4.1U ) gener- 



ates precisely such a velocity. The density p(x) = p(x, t = 1) that we shall find in the course 
of this proof will then be the matrix chain eigenvalue density at the Kosterlitz-Thouless 
point. 

The velocity condition (|4.11|) can be checked conveniently using the implicit Hopf 



solution ( |2.27| ). We start at t = 1/2 with the Hopf function 

XxJ 



fi/2(x) = ^\og(J-)+0(x 3 ) (4.12) 
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and we must finish at t = 1 with / being equal 



= f(x, t = l) = 2x + mp{x) + 0{x c 



(4.13) 



Since the elapsing time interval is At = 1/2 the implicit solution ( |2.27| ) imposes the 
following relationship between fi/ 2 and /i 



fi(x) = fi/n[x - \fi(x)]. 



(4.14) 



This is a complex valued equation which has both real and imaginary parts. However, 
the only unknown in fi/ 2 and f± is one real valued function p(x). Therefore, the internal 
consistency of equation (|4.14| ) is not automatically guaranteed. Rather, such consistency 
is an indication of the correct choice for pi/2(x) that we made in ( f4.10| ). 

For the purpose of computational convenience let us represent p(x) in the form 



Trp(x) = ir\x\a(x) + 0(x 3 ). 



(4.15) 



Note that again the x < and x > regions do not interact. Therefore, we can safely 
replace \x\ — ■> x in ( |4.15| ) so long as we apply our answers only to x > 0. Then taking 
( [4.12|) and ( [4. 13| ) into account, equation ( fll4| ) translates into 



2x + iiTxa(x) = —2xa(x) log 



— — Xxa(x) 



Using log(— iy) = logy — in/2 we see that the imaginary part of this equation is satisfied 
identically. The real part yields then an equation on cr(x), 



2x = —2xa(x) log 



— Xxa(x) 



or, equivalent ly, 



— *— + log a(x) = log( ^- ) 
a{x) \Ax J 



(4.16) 



with A = n\/2. For small x the solution of ( fll6| ) is approximately 

cr(x) oc 



log [1/(A 



x 



in complete agreement with ( |4.1|) . 

A remarkable property of equation (|4.16|) is its universality. Indeed, the interaction 
potential enters p(x) only through the 0(x 3 ) terms which did not matter in our analysis. 



Therefore, the logarithmic terms summarized by cr(x) will remain the same for all, not only 
quartic, matrix interactions. Such universality may indicate that equation ( 4.1(j| ) contains 



certain information about the continuum limit of the theory. To understand what precisely 
this information is would be very interesting. A more specific problem that can hopefully 
be addressed using ( f4.16|) is to find the free energy at e = 1. Once we know the critical 
index of the free energy 7 str we can find the effective central charge at e = 1 and determine 
what kind of a conformal field theory the matrix chain describes at that point. If the 
relation between 7 str and the eigenvalue index 5, 

7s tr = 1-6 

holds in this case, formula indicates that at e = 1 the string susceptibility is most 
likely zero. This is characteristic of c = 1 theories. However, the logarithmic corrections 
to the free energy at e = 1 would certainly be different from the logarithmic corrections in 
the c = 1 theory on a straight continuous line. Therefore, these corrections deserve more 
investigation. For example, equation (|4.16|) may have an interpretation in terms of vortices 



which start forming a strongly interacting system precisely at the point e = 1. 

To compute the free energy carefully one must know the eigenvalue density in a certain 
interval of coupling deviations 5g = g — g cr . For small values of 5g the eigenvalue density 
can be found, as before, using the method of dispersion relations. 

Again, we shall look for p(x) which has two close square root branch points at x ~ 
±i\/b and is an even real function for real values of x. We shall seek an approximation to 
p valid for x ~ Vb <C a where a is the typical width of the quartic potential. Moreover, 
since the density Pi/2(x) seems to have a simpler form than p(x) itself we shall continue 
working directly with pi/2(x) even for 6^0. 

The eigenvalue density with the desired properties can be found as a series in powers 
of rj = y/b/x. That is, we go to the region \/b <x<a which is still within the limits of 
our approximation. In this region r\ <C 1 is a good expansion parameter. Consequently, we 
may look for p\i%[x) in terms of a series in powers of r/ or, equivalently, of b. We shall now 
demonstrate that the relevant, consistent with the hydrodynamic picture series is given by 

Pi/2(x)=p 1/2 (x)\ b=0 + J2 a kT^ k -^- (4-17) 

k=l 1 1 

This ansatz for pi/2(x) is not arbitrary. A similar expansion arises, for example, in the 
e < 1 phase of the matrix chain. There p(x) ~ \/b + a 2 x 2 can also be expanded in a series 

I ^ 00 /l/2\ b k 

Mx)=a\x\^l + ^- 2 =a\x\ + ^ k 



\2k-l 



identical in form to ([4 . 1 7| ) . Both of these expansions are explicitly singular at x = but 
the singularity disappears when the full series is summed. 

Most interestingly, expansion ( [4. 171 ) is consistent with the Hopf equation and the 
boundary conditions for any set of real coefficients a&. To see this we simply repeat the 
computation of equations ( |4.12|) — (|4 . 1 6|) . The function fi(x) = f(x,t = 1) is still given by 
equation ( [4. 13| ) while for fi/2(%) we now have 

+ in > n jt -^r-_ j + 0(x 3 ). 
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(4.18) 



Imposing the implicit relation ( J4.14 ) we find that again the imaginary part of ( 4.14| ) is 
satisfied identically while the real part yields an equation on p(x), 



2x = — 2p(x) log 



~\p(x) 
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(4.19) 



The coefficients a& have been left undetermined by the Hopf equation. Rather, they are 
constrained by the analytic structure of /Oi/ 2 (x). Indeed, if pi/2{x) has two cuts shown in 
fig. 3 then the real part of pi/2 on the cut edges follows directly from ( |4.17| ). For example, 
on the right edge of the upper cut, where x = iu + 



r{iu 



= Re [7rpi/ 2 (x)] = Re 



— log! 

7T V AX 



(4.20) 



2u. 



Using the square root nature of the branch points we can then perform an analytic con- 
tinuation and find the real part of pi/ 2 on the other three edges. 

The whole 7rpi/ 2 (^) can now be restored through a subtracted dispersion relation as 
in (|3.40|) . The contour of integration is the same contour C in fig. 3 while the dispersion 
integral itself is a little different, 

dZ (4.21) 
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The integrand in ( fl.21|) has been constructed in such a way that the contribution from 
infinite loops closing the contour C vanishes (see ( [4.17]) ) while the cuts contribute only 
through the real part of 7rpi/ 2 (^), 



2i r(iu) 
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On the other hand, the integral in ( f4.21|) can be evaluated via residues to give 
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where r = ^P\ii{x — 0). Comparing this to Q4.22| ) and evaluating the integrals we find 
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(4.24) 



The constant ro, similarly to in (|3.44|) cannot be determined by considering the region 
x ~ Vb < a. To fix r one would need to construct pi/ 2 for all, not only small values 
of x. The analytic structure of such globally constructed p(x) would include two other 
branch points at x = ±xq ~ a. Indeed, in matrix models with polynomial potentials the 
eigenvalue density is usually nonzero only within a finite interval [— xq, xq]. The endpoints 
xq are, as a rule, square root branch points. From the viewpoint of small x expansions, the 
square root character of the outer branch points is a nontrivial property of the eigenvalue 
distribution which is not automatically guaranteed. Therefore, imposing it is likely to 
constrain the additional parameters such as ro in Ql24| ) or r 2 in ( |3.44j ). After all, a very 
similar phenomenon has occured when we imposed the two cut structure on pi/ 2 {x). There 
requiring two square root branch points at x = ±iVb fixed at once an infinite number of 
coefficients which remained unconstrained by the Hopf equation. 

To summarize, when the Hopf equation or the functional equation G_|_[G_(x)] = x 
are combined with reasonable assumptions about the analytic structure of the eigenvalue 
density the solution becomes determined uniquely. In the next section we shall explain, 
using an example, how such global solutions can be constructed. 



5. Exactly Solvable Interaction Potentials 

By now we have seen that the Euler equations, the Hopf equation and the functional 
equation ( |2.29|) are all very useful in studying the singularities of p(x). In this section we 
shall show how the same equations can be utilized to construct certain "exactly solvable" 
potentials V(M) for which p{x) is an elementary function. We shall demonstrate that for 
these potentials the expansion methods of sections 3 and 4 do reproduce the exact results. 

The basic idea of such exact solutions is very simple. One chooses at will a midway 
eigenvalue density ^1/2(2^) making sure it is properly normalized. Then the fluid is allowed 
to evolve under Euler equations. Finally, at t = 1 one reads off the velocity v (x, t = 1) and 



the corresponding eigenvalue density p(x,t = 1) = p(x). The result is an exact solution — 
given by p(x) — for the model where the potential can be found from (|3~T|) , ( pTTOj ) 



U'(x) =2[x-v(x,t= 1)]. (5.1) 

The hardest part of the problem is to choose p\/2(x) correctly, so that the resulting U (x) is 
physically reasonable. Here we shall present a construction ]16| which produces a Penner- 
type double-well potential rather similar to the pure quartic potential of sections 3 and 4. 
To this effect consider a special pi/2(x) defined by the following quadratic equation 

k(x 2 + r 2 ) 2 + m 2 (x 2 + r 2 ) - 2(x 2 - r 2 ) + b = (5.2) 



where r = 7rpi/2(x)/2. This equation looks, of course, very artificial, and it is. As it 
turns out, the potential U(x) corresponding to such a choice of pi/ 2 is particularly simple. 
Later on, we shall replace the left hand side of (|5.2| ) by a more general polynomial to get 
a solution for the e > 1 phase of the matrix chain. 



The potential U and the eigenvalue density p that follow from ( |5.2| ) can be easily 
found. As in ( |4 . 1 2| ) — ( 4.14|) we use the implicit Hopf solution 



h(x) = A/2 [a: - \fi{x)] 



(5.3) 



with fi/2(x) = inpi/ 2 {x) and 



A(^) = x U'(x) + inp(x). 



(5.4) 



To determine A we replace x — > x' = x — f\(x)/2 and r(x) — > r(x') = 7r / o 1 / 2 (x / )/2 in 
equation 05.2ft . As a consequence of (|5.3[) 



r(x') 



which we can substitute into ( |5.2|) to get a quadratic equation on A 



(1 - kx 2 )(x- A) 2 -m 2 x(x- A) -b + x 2 = 0. 



(5.5) 



Solving this quadratic equation we obtain the final value of the function / 



and, therefore 

m 2 x 




1 kx " (5.7) 

' ' +(b-x 2 )(l-kx 2 ). 
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The original definition of the matrix chain deals, however, with a rescaled potential V(M) 
related to U(x) by (|2.9| ). In our model 

where we introduced the "scaled" mass k and the coupling constant G 

(5.9) 



m 2 = 2 + k 2 6 2 



k = ±Ge 3 . 



As always, by rescaling all matrices in the matrix chain partition function ( |1 . If ) we can 
make k 2 = V"(0) anything we want. Changing the value of k 2 simply changes the scale 
of e and shifts the position of all critical points. We shall adhere to the convention of 
( p.31|) that was used for pure quartic potentials and set k 2 = —4. In terms of the original 



parameters m 2 and k this would mean 

! = -2(2e 2 - 1) 

(5.10) 



k = ±Ge 3 , 

a relation quite similar to ( |2.32|) . 

The resulting potential V(x) is plotted in fig. 4. It goes to infinity as x — > ±x r 
^/2/e 2 G thereby restricting the eigenvalues to a finite range [— x max , x max ]. Note that the 
plot of fig. 4 makes sense only when K 2 x max + Gx max > — otherwise for x — ► ±x max the 
potential goes to minus, not plus infinity and therefore is unbounded from below. For 



j max 



K 



2 



—4 this constraint translates into e < l/y/2 or, equivalently, m 2 > 0. 



Such a picture is certainly consistent with what we already know. Indeed, at e = l/v^2 
we expect p(x) to develop a logarithmic singularity. On the other hand, the solution (|5.7Q 
cannot, by construction, have such logarithms. Therefore, it is perfectly natural that this 
solution should degenerate precisely at e = 1/ a/2. Furthermore, our exact solution has no 
singularities for e between and 1 / v2. Here is an example of a special situation when the 
terms responsible for such singularities — say, the terms in the curly brackets in (|3.25|) — 
are arranged to vanish. In fact, the exact eigenvalue density (|5~77l ) can be used to verify and 
confirm all of the perturbative techniques developed in section 3. For instance, when 6 = 




adapting the small x asymptotic series of ( |3.19| ) to the case of our nonpolynomial potential 
yields precisely the expansion of the exact solution ( |5.7|) in powers of x. Moreover, for 
b 7^ it is possible to calculate, up to arbitrarily high orders, the small b expansion similar 
to ( |3.3U| ). The result, not surprisingly, can be summed into an exact formula thereby 
providing a nontrivial check of our computational methods. 



Quite remarkably, the same exact solution (|5.7f) can be used to rediscover not only 
the singularity at e = 1/ \[2 but also the Kosterlitz-Thouless phase transition. To see this 
consider the same potential ( |5.8|) but with k 2 > and G < 0. Such a potential, also 
sketched in fig. 4, has two humps at Xh = ±k/ \/\G\ which play exactly the same role as 
the single hump at x = 0. By analogy with ( |2.31| ) we shall fix the second derivative of the 
potential at the hump to be V"{x-^) = —4. This imposes the following restriction on k 

k 2 = -. (5.11) 

1 — t z 

That is, the potential degenerates as we approach the Kosterlitz-Thouless phase transition 
at e = 1. At that point expression fl5.7Q cannot possibly give the correct eigenvalue density. 
Such "nontranscendental" form of p does not contain the requisite logarithmic singularity 
which, as we found, has to be present for any, even nonpolynomial, interaction potential. 

Most likely, one could use the same methods to construct the exact solutions that 
would describe the logarithmic singularities themselves. The idea would be to invent a 



globally defined Pi/ 2 {x) having both the small x and small b expansions consistent, to first 
order, with (|3.23|) , ([3.44] ) and ( 4.24 ). Then evolving it to t = 1 would, by construction, 
generate the p(x) with an adequate amount of logarithms. After the global p(x) has 
been found the direct computation of the free energy and various critical indices becomes 
reasonably straightforward. It would certainly be most interesting to actually carry out 
this program, especially for the Kosterlitz-Thouless phase transition where the results 
would have a direct meaning in terms of string theory. 



6. Interaction of Vortices with Curvature Defects 

We have seen throughout this paper that the Kosterlitz-Thouless phase transition and 
the subleading logarithmic singularities of the e < 1 phase bear remarkable mathematical 
similarity. This similarity is perhaps expressed best by the formulas for pi/ 2 (x). Indeed, 
an explicit expression for p 1 / 2 at e = l/y/2, 

Tvp 1/2 (x) = a/26 + Ax 2 + gu(x) (6.1) 



with 



t \ 8 U 2 b \ / 9 & , 

u(x) =—x \2x ^ 2 J \ X 2 



-2r 2 (2x 2 + ^\ -b(2x 2 + ^ 



\/x 2 + b/2 + x 
x 2 + b/2 — x 



(6.2) 



has precisely the same structure as its analogue for e = 1, equation (|4.24| ). 

Such parallels suggest that the qualitative physical picture of logarithmic singularities 
must be related to the vortex picture of the Kosterlitz-Thouless transition. Below we shall 
argue that this is indeed the case. 

In continuum theory, the action for a one-dimensional noncritical string compactified 
on a circle of radius R is 

S[g,X] = ^J d 2 x^gg^d a Xd p X (6.3) 

where X is a compactified string coordinate, x = X + 2%nR for any n G 2 and g a p - 
the two-dimensional worldsheet metric. For flat g a p a typical vortex — a solution of the 
equations of motion with a nonzero winding number n — is given by 



X(r, 4>) = n(j>R 



(6.4) 



r and <p being the polar coordinates parametrizing the flat area of the worldsheet. The 
value of the action for such a vortex would be infinite were it not for the cutoffs. Let 
us therefore introduce the infrared cutoff L and the ultraviolet cutoff a (say, the lattice 
spacing if the worldsheet is discretized.) With these cutoffs the vortex action equals 

S vort = 2n 2 R 2 \og[^j. (6.5) 

The total weight with which a vortex contributes to the partition function of the whole 
theory can be found by multiplying the number of places where the vortex could be centered 
(which is simply the total number of lattice sites iVtot ~ (L/a) 2 ) by the Gibbs factor 
exp(- SV^) 

{ L\ 2 / l n 2(l-n 2 fl 2 ) 

w n oc - e _s ™ rt = (-] ■ (6-6) 



As a result, depending on the value of R vortices may or may not be important |]TD] . If 
R > 1 one has 2(1 — n 2 R 2 ) < for all integer n and the total weight of any vortex vanishes 
when the cutoffs are removed. In this case vortices can be neglected. On the contrary, for 
R < 1 the weight wi goes to infinity whenever L/a — > oo and so vortices do matter. The 
two phases which arise are separated by a Kosterlitz-Thouless phase transition with the 
critical value of compactification radius R CT = 1. 

Let us now imagine that the two-dimensional worldsheet is not completely flat but can 
have freely moving conical singularities. Around a conical singularity such a worldsheet 
can still be described by the polar coordinates r and <p but with <p running in a different 
interval < <fi < 2n + a. The "excess angle" a characterizes the amount of curvature 
acquired by the manifold due to the conical singularity. 

For a vortex centered right at the top of the cone the vortex field X would be given 

by 

X{r,<j>) = -^ncj>R. (6.7) 
2-7T + a 

Consequently, the vortex action changes as well, 

Scone = -^—n 2 R 2 log (-). (6.8) 
2-7T + a \ a J 

As we see, if a > — that is, when the conical singularity contributes negative curvature — 
the value of the vortex action decreases compared to flat space. Furthermore, in two- 
dimensional gravity it does not cost any extra energy to create a conical singularity, so 
long as the total worldsheet area and the genus remain the same. Therefore, we could 
contemplate a Kosterlitz-Thouless phase transition induced by such pairs of "correlated" 



cones and vortices centered closely to each other. The critical value of R or the corre- 
sponding e = 1/R for this transition can be found similarly to ( |6.6| ), 



V Z7T + a 

Of course, the pairs of cones and vortices with nearby centers form only a small susbset 
of all possible states in our system!. Roughly, the numbers of possible positions of either 
a vortex or a conical tip are (L/a) 2 each, adding up to a total of ~ (L/a) 4 states. Only 
~ (L/a) 2 of them are the states where the cone and the vortex have the same center. 
Consequently, the contribution of such "correlated pairs" can produce at most a finite size 
correction to the free energy which is down by (a/L) 2 oc 1/A for a worldsheet of area A. 
But this is exactly what we find from (|3.23| ). There the logarithmic correction to p(x) 
is also subleading. Indeed, at e = l/y/2 the leading contribution to the free energy — the 
one which is important in the double scaling limit — is associated with the first term of 

expansion for the eigenvalue density, p(x) cx \x\ + For any e < 1 this contribution is 

simply the standard c = 1 free energy, 

F(5g)K-^— + 0(5g 3 ). (6.10) 
I logdc/| 

The logarithmic correction x 3 log[l/(Ax)] would modify only the 0(5g 3 ) terms of ( |6.10| ). 
And, as it turns out, such 0(5g 3 ) terms correspond precisely to 1/A corrections. 

To see this, we just have to recall that the cosmological constant 5g and the worldsheet 
area A are a pair of thermodynamically conjugate variables. That is, the string theory 
partition function with a nonzero 8g, given by J r (5g), is related to the partition function 
of random surfaces with fixed total area Z s t r (A) by the formula 



Z str (Sg) = F(6g) = J dAe- S9A Z str (A)- 



Consequently, an expansion in powers of 8g arising in (|6.10|) would translate into a 1/A 



expansion for the fixed area partition function and the fixed area free energy. 



5 Here is, perhaps, the most objectionable part of this argument. Strictly speaking, the 
Kosterlitz-Thouless picture can be used reliably to infer the existence of a phase transition only 
for the whole system, and not for an artificially chosen small subset of states. The only justifi- 
cation for what we say is that the results gotten via such reasoning do not contradict the exact 
computations of sections 3 and 4. 



Finally, this picture can be used to explain why the logarithmic singularities form a 
discrete set. The reason is, for a discretized surface generated via a matrix model the 
excess angle a at any graph vertex can assume only certain discrete values. For example, 
a matrix model with a cubic potential would produce surfaces composed of equilateral 
triangles. On such surface the angle a corresponding to a vertex with the coordination 
number q (where q triangles meet together) equals 

« = 3(9- 6). 



We see that a and therefore the critical lattice spacings computed from (|6.9|) do come out 
discrete. However, this rough estimate fails to yield the exact critical values of e. This 
should not be a surprise. To get these critical values right one would have to take into 
account various fluctuations that we neglected. For instance, the picture of an almost flat 
lattice with just a few curvature defects scattered around is certainly quite inaccurate. In 
reality, almost every vertex of a generic two-dimensional graph would have q ^ 6. It would 
be quite interesting to find out whether such effects can be taken into account and whether 
our rough picture survives that. 

To summarize, the logarithmic singularities at e < 1 carry along lattice information. 
That is, their position and order depend on the particular type of polygons tiling our 
random surface. However, this does not yet mean that such effects are of little interest. 
For example, the Kazakov's multicritical points in one-matrix models [^2j provide, also 
through a very special way of tiling, a description of minimal models coupled to quantum 
gravity. There the connection with continuum theory is not obvious from the matrix model 
formulation. A similar indirect continuum interpretation might exist for the logarithmic 
singularities in the matrix chain. In any case, the subleading singularities of the e < 1 
phase are quite an unusual phenomenon for matrix models. Indeed, in most theories the 
whole singularity structure is captured by the double scaling limit of the matrix model, the 
corrections to this limit being perfectly regular. For the matrix chain this does not seem 
to be so, which perhaps is rather satisfactory. As an example, it is very well known [[F] that 
the double scaling limit of the chain model is smooth throughout the interval < e < 1 
and describes c = 1 string theory. This leads to a surprising conclusion that the discrete 
structure of the target space has no effect at all. We see that although such a conclusion is 
absolutely true in the continuum theory, the target space discreteness does show through 
a nontrivial set of corrections. 

Finally, the logarithmic singularities can be of interest for a separate reason. It appears 
remarkable that such a complicated family of singularities arises precisely on the boundary 
between the c < 1 and c > 1 theories. One might wonder, maybe somewhat groundlessly, 
whether a similar kind of structure could also emerge in any of the possible c > 1 models. 



7. Conclusions 

We have investigated, in the large iV limit, the phase structure and the eigenvalue 
density of the one-dimensional infinite random matrix chain. We developed a systematic 
expansion applicable in the vicinity of the Kosterlitz-Thouless phase transition as well as 
throughout the e < 1 phase of the model. Remarkably, such expansion can be used to 
extract the critical properties of the eigenvalue density even when the model is hard to 
solve exactly. 

It turned out that the infinite matrix chain exhibits a few effects rather unusual in 
matrix models. In particular, we found an infinite number of subleading (not important 
in the double scaling limit) critical points, similar in nature to the Kosterlitz-Thouless 
phase transition. These special points appear to arise from interactions of vortices with 
curvature defects and depend on the specific type of polygons tiling a random surface. 

The computational tools that were required here, such as the hydrodynamic rep- 
resentation and the functional equation G+[G-(x)] = x are likely to have many more 
applications. Even in the case of the matrix chain these methods reproduce and generalize 
with great ease the results which would look nontrivial from other viewpoints. They can 
also be used to explore the e > 1 phase and to establish the parallels with the character 
expansions. 

In addition, there are a number of specific open problems that should be possible 
to resolve with our tools. One such problem is to compute carefully the free energy of 
the matrix chain at the Kosterlitz-Thouless point. This could be done via the procedure 
outlined at the end of section 5. Furthermore, it must be possible to relate the critical 
behavior of the eigenvalue density p(x) that we now know, to the physical picture of vortex 
interactions. Finally, it would be very interesting to evaluate, or at least to represent in 
terms of the Hopf equation, the correlation functions of the matrix chain. Such repre- 
sentation would help to separate and compute the contributions of vortices with different 
vortex charges. As a consequence of this computation, we shall be able to study in a very 
direct fashion the dynamics of vortices interacting with quantum gravity. 
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